Optimized Double-well quantum interferometry with Gaussian-squeezed states 
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A Mach-Zender interferometer with a gaussian number-difference squeezed input state can exhibit 
sub-shot-noise phase resolution over a large phase-interval. We obtain the optimal level of squeezing 
for a given phase- interval A^o and particle number N ^ with the resulting phase-estimation uncer- 
tainty smoothly approaching 3.b/N as A^o approaches lO/A/", achieved with highly squeezed states 
near the Fock regime. We then analyze an adaptive measurement scheme which allows any phase 
on (— 7r/2,7r/2) to be measured with a precision of 3.5/A^ requiring only a few measurements, even 
for very large N . We obtain an asymptotic scaling law of A^ (2.1 + 3.2 ln(ln(A^tot tan A^o)))/A^tot, 
resulting in a final precision of ^ 10 /Ntot- This scheme can be readily implemented in a double- 
well Bose-Einstein condensate system, as the optimal input states can be obtained by adiabatic 
manipulation of the double- well ground state. 

PACS numbers: 03.75.Dg,03.75.Lm,42.50.Dv 



Measuring an arbitrary phase with precision well above 
the standard quantum limit (SQL), i.e., below shot-noise, 
has been a long-standing challenge in quantum inter- 
ferometry [il, i, S, Q. The SQL minimum phase un- 
certainty is 1/y/Ntot^ but the theoretical lower limit to 
the phase uncertainty, known as the Heisenberg limit 
(HL), is 1/Ntot^ where Ntot is the total number of par- 
ticles used in the determination of the phase. There 
have been many proposals to achieve I /Ntot scaling in 
a two-input interferometer, which are based on number- 
difference squeezed input states [l, Jsl, [HI, coherent 
and/or squeezed vacuum input states p, [3, |8|, or the 
maximally- entangled A/'-particle NOON state, ^{\N^0) -\- 
\0,N)) [1 E El [12. Recently, the double- well Bose- 
Einstein condensate (BEC) has emerged as a promis- 
ing system for hi gh-p recision matter- wave interferome- 
try l3, 14, 15, 1^, 17], including progress towards atom- 
counting at the single-particle level [18, 19]. For this 
system, the squeezed-vacuum protocols are not appli- 
cable, while the NOON state is not suited to deter- 
mine an unknown phase due to the periodicity of the 
phase-distribution 0,El]- This leaves number-difference 
squeezed states as a viable candidate, although to date 
there has been no systematic study of how to measure 
arbitrary phases at or near HL precision in the large- 
limit. In this Letter we perform such an analysis 
and show that an asymptotic scaling of {\n{\n Ntot)) / Ntot 
can be achieved via multiple adaptive measurements with 
Gaussian number-squeezed states, which can be readily 
created in a double- well BEC. 

For measuring a phase of ^ = 0, it has been shown that 
the Twin-Fock (TF) state, and the related Pezze-Smerzi 
(PS) state can achieve Heisenberg scaling [20',^, 23, 23 |. 
We find, however, for ^ 7^ the phase-uncertainty of the 
TF and PS states rapidly decay to worse-than-SQL, and 
in the limit of large N approach constant values, inde- 
pendent of A^. While the PS state was only investigated 
for l9 = 0, Kim et al [24] investigated l9 7^ for the TF 
state with N = 100. They claim a phase uncertainty 



~ I /Ntot for < l/N^ and growing rapidly thereafter. 
Our results similarly indicate that the TF and PS states 
become worse than shot-noise for ^ 1/N. 

To find the optimal input state, we constrain ourselves 
to the ground-states of a double- well BEC with repulsive 
interactions for experimental obtainability. The double- 
well BEC system is described by the Hamiltonian, 
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where r is the inter-well tunneling rate, U is the atom- 
atom interaction strength, and S is the asymmetric tilt of 
the double-well, presumably due to the external pertur- 
bation being measured. The angular momentum opera- 
tors are defined as Jx = ^{c\^cr -\-c\^cl), Jy = ^{c\^cr — 
c\^cl) and Jz = ^{c]^cl — c\^cr), with cl^cr being the 
annihilation operators for particles in the two localized 
modes. For repulsive atom-atom interaction, and J ~ 
the ground state is very close to a Gaussian squeezed 
(GS) state of the form \a) ex Eni-Ar/2 ^"'''^'^''V). 
where \n) is a number-difference eigenstate satisfying 
Jz\n) = n\n). The width a, depends on the param- 
eter u U/r, and is given by 7V/4V1 + uN |25|. 
The nature of our adaptive measurement scheme requires 
that we tune u to an optimal value which takes into ac- 
count our prior knowledge of N and ^, thus we have 
u u{N^O). This tuning is accomplished by varying 
U/r via a Feshbach resonance and/or changing the shape 
of the double- well potential, and allows a to be varied be- 
tween and a/7V/2, corresponding to maximal number- 
difference squeezing and no squeezing, respectively. 

To implement a Mach-Zehnder interferometer (MZI), 
we set U = and allow tunneling for t = 7r/4r duration, 
thus realizing a linear 50/50 beamsplitter, described by 
the propagator e*^'^^. This is followed by a sudden rais- 
ing of the potential barrier to turn off tunneling and allow 
phase acquisition due to the small but non- vanishing S. 
Holding the system for a measurement time T, a phase 
shift of ^ = —ST will be acquired, described by the prop- 
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agator e*^^^. The barrier is then lowered again to imple- 
ment a second beamsplitter. 

In the MZI, a symmetric input state |^in) is trans- 
formed into the phase-dependent output state |^ont) = 
gif j.gi^j.gif j^l^.^^^ = e-^^^^e^^^-|^in) [1]. Apply- 
ing this transformation to a typical GS state, results 
in an output state whose properties are easily under- 
stood using the Bloch sphere quasi-probability distri- 
bution 0, which assigns a probability to each point 
on a sphere of radius N/2 according to P{0^(j)) = 

|(Ar/2|e^'^^(^/2-^)e^^^-^|^ont)p. Note any TV-atoms state 
for which all atoms are in the same single-particle state 
can be written in the form e-^^-^e-^^^(^/2-^)|7V/2), 
which has eigenvalue N/2 with respect to the projection 
of J onto the axis defined by O^cj). In this picture, the 
TF and PS states are thin equatorial rings, and a GS 
input state is an ellipse centered on the Jx axis, com- 
pressed along the J^-direction. Typical GS input and 
output states are shown in figure 1(a). 

Phase information is obtained by measuring the num- 
ber difference between the two interferometer modes, 
which projects the output state onto a eigenstate. 
Quantum fiuctuations in this measurement are governed 
by the projection of the output distribution onto the J^- 
axis. Due to the rigid rotation, the width of the projec- 
tion will be determined by a ^-dependent combination 
of the Jz and Jx noise of the input distribution. The 
goal of this paper is thus to find the optimal amount of 
squeezing to minimize the phase uncertainty given a fixed 
particle number N and an initial estimated phase with 
uncertainty A^o- 

Before we present numerical results from a rigorous 
Bayesian analysis, we first use linearized error propa- 
gation to provide an approximate analytical description 
of the interferometer performance. An analytical re- 
sult is important for predicting the behavior at large 
N^ where a numerical result is inaccessible. In this ap- 
proach, the phase uncertainty is estimated by evaluating 
A6 = [d{Jz) /d6]~^AJz at the interferometer output. For 
input states symmetric around n = 0, the expectation 
values at the output are related to those at the input via 

{I 



sin l9(4) and AJz = \J cos^ OAjf ^ sm^ OAjf . 

For GS states, (J^) ^ N/2, and AJ^ = cr, which im- 
mediately leads to {Jz) = N{sm6)/2. In figure 1(c) we 

N/2 



see that AJi 



Since GS state 



is a minimum uncertainty state with AJ*AJ* = (J*)/2, 
we see that AJ^ = N/4a. This leads to AJ^ = aN/a^, 
with a ~ 0.06. Exact numerical calculations verify this 
analytic form for 1 <C cr <C VN/2, as shown in figure 
1(d), but with a = 0.09. Inserting these results into the 
error-propagation formula, we find 




0.09A^ tan l9 




FIG. 1: (Color online) Bloch-sphere analysis of the MZI with 
a GS input state: (a) A typical GS state at input and output 
stages; (b) output states for optimized input states with = 
0,±Y^,±^,zb^,±^,±-||; (c) geometric origin of the Jx input 
noise, AJ*; (d) numerical results plotting AJl/N versus a 
on a log- log scale for three different A^- values, validating the 
functional form of A J* derived geometrically from (c). The 
dashed vertical lines correspond to a = where A J* drops 
to zero. 



The TF and PS states roughly correspond to a fixed a ^ 
1, resulting in limAr^oo A^ = .18tan^/<j^, which quickly 
becomes saturated to an A/'-independent constant for 9 ^ 
1/A^, a result we have verified numerically with exact 
Bayesian calculations. On the other hand, if holding u 
fixed so that a ~ A^^/^, the phase uncertainty scales as 
A^ ~ l/A/"^/^ for ^ = 0, as discussed in [26], and would 
eventually saturate to ~ I/a/TV for ^ 7^ 0. Rather than 
holding (J 01 u fixed, we propose varying u and thus a 
with N in order to minimize the phase variance. By 
setting dAO/da = we find 



A^, 



r{0,N) 
^{0,N) 



.503(7Vtan|l9|)^/^ 
1.23(tan|^|)^/V^'^^ 



(3) 

(4) 



(2) 



From self-consistency, these expressions are valid only 
when 10/N < |6>| < tan-i(.137ViV) ^ 7r/2. 

We now employ rigorous Bayesian analysis to quan- 
tify the phase uncertainty and validate our approximate 
analytic results, again assuming that \0\ is not too close 
to 7r/2. According to Bayes theorem, upon a measure- 
ment result nm, the probability that the actual phase is 
is P((^|n^) = P(n^|0)/ / deP{nm\0), where P{n\e) = 
|(n|e~*^^^ The error-propagation result (|4]) is very 

close to the 68% confidence interval of P(^|n) because 
the underlying number distribution of the optimized GS 
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output state is well approximated by the gaussian distri- 
bution P{n\0) 

An = 



1 



2tan2 6>c 



[V^An]-ie-("-^^i"^/2)/2A- , where 
^minifiaiN)^ with being the un- 



tan2 e 



known phase, and Oa being the assumed phase used for 
optimization. Provided that \0\ ~ |^a|, the dependence 
on 9 is weak, and we have An ^ ^/^j2Gmin{0a-, N). 
This shows that the most-probable outcome is n = 
A/'sin^/2, which is sensitive to the sign of 0. Because 
An is only weakly dependent on ^, the inverted dis- 
tribution P{(j)\n) will also be close to Gaussian in the 
small-angle regime, with width given by (|4|). To make 
a theoretical performance analysis for a fixed ^, we av- 
erage over all possible measurement outcomes, defining 

^(^1^) = En^(^l^)^(^l<^) E 0- This can be in- 
terpreted as the probability for an experimenter to infer 
(j) given a true phase-shift of 0. The phase uncertainty 
A^ is then defined as the 68% confidence interval, via 

Using this approach, together with the exact double- 
well ground-state, we numerically find Umin^ the value 
of U/r which minimizes the phase uncertainty. In figure 
2(a), we plot the corresponding (imin = o-{umin) as a 
function of N for several ^s. Also shown is a least-squares 
fit to the N > 10^ data (including many data points not 
shown explicitly) to the analytic form (|3]), giving 



^ ( 1.00, 
\ 0.45( 



1^1 < 10/iV; 
45(7Vtan|(9|)l/^ \0\ > 10/N, 



(5) 



in good agreement with our analytical result. Inverting 
Eq. ([5j) to leads to 



N 



52 1_ 

{tan \e\y/3 1^1/3 Ar> 



1^1 < 10/iV; 
1^1 > 



(6) 



In Fig. El^b) we plot the corresponding minimized 
^^min versus N for several phases, achieved by setting 
o' = (Jmini^j N)' Again fitting the N > 10^ data to the 
analytic form of (|4]), we find 



r 3.50/N, \e\ < 10/N; 

M.63(tan|e|)i/3/iv2/3, \0\>1O/N. 



(7) 



The difference between the prefactor here and (|4]) is pri- 
marily due to a factor of approximately \/2 which comes 
from the definition of P{(j)\0). 

In practice, is not known a-priori, hence it is not 
clear what value for to use in determining amin{O^N) 
via Eq. ([5j). If we assume prior knowledge of the 
form P{0) (X exp[-(l9 - 6>o)V2A(9g], we should first re- 
move ^0 by adding Oq/T to the tilt 5 during phase ac- 
quisition, and then use amini^Oo, N). After obtain- 
ing a measurement result ni, the estimated uncertainty 
A^i is then be determined via Jq^^^q^ dO' P{6'\ni)^ with 
P{0'\ni) being given by Bayes theorem. This will result 
in A6>i ^ A0min{^0o,N). Based on Eq. (O, this un- 
certainty appears to scale only as A^~^/^, only a slight 
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FIG. 2: Figure (a): Optimal width amin versus N for different 
0; (b): corresponding minimized phase uncertainty AOmin- 
In both figures, from (i) to (iv) the interferometer phases are 
^ = 0, 0.01, 0.1 and 1. The data points represent numerical 
results from strict Bayesian analysis using the exact ground 
states of a double- well BEG, while the straight lines represent 
the asymptotical forms of Eq.(|5| and Eq.(|7|, respectively. 



improvement of N'^^^ over the SQL. However, in many 
applications requiring high precision, the phases are very 
small, in which case the phase uncertainty can be reduced 
considerably due to the explicit phase dependence in ([7j). 
This is in contrast to a shot-noise-limited interferome- 
ter, where A^ = for all ^s not too close to ±7r/2. 
The explicit theta-dependence in the optimized scheme 
is due to the fact that stronger number-squeezing can be 
tolerated at smaller angles before the J* noise becomes 
detrimental. For example, if the phase is known to be 
smaller than 1/V^, we have AOmin < 1.63/7V^/^, which 
is now an N'^^^ improvement over the SQL. As can be 
seen from Eq. ([7|), a maximum sensitivity of 3.5 /N can 
be achieved for |^| < which is true Heisenberg 
scaling. 

In fact, almost any phase between — 7r/2 and 7r/2 
can be measured at the maximum precision of 3.5 /N 
if the present scheme is combined with multiple adap- 
tive measurements 'tJ. After the first measurement 
as described above, we can again rebalance the interfer- 
ometer by adding 61 /T to the tilt, followed by a sec- 
ond measurement with (72 = amini^^i^ with result 
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722. The Bayesian distribution for O2 will then be ap- 
proximately P(6>2|^2,^i) 0^ exp[— (^2 — '^2)^/2^2] 7 where 
-^(2nj/N) and I/t/^ - 1 /a^2 



sm 



ELol/A^i. Since A^, 
is much smaller than AOj-i we can say yj ~ A6>j ~ 
/\Omin{^^j-iiN). In other words, since the distribu- 
tion after a measurement is much narrower than the 
previous distribution, multiplying the distributions has 
little effect, so that the final uncertainty is effectively 
determined by the resolution of the final measurement 
alone. After M iterations, with Gj = cFmin{^^j-i-,N) ^ 

0.57(7VtanAi9o/2.1)^"', we find 

Mm - (2.1/^") (TVtan Ai9o/2.1)^"'^ . (8) 

While the above expressions are good estimates of the 
expected behavior, in practice each cFj and A6>j would 
be computed exactly by applying Bayes theorem after 
each measurement. This procedure should be repeated 
only until A6>m ^ 10/A^, after which an addition mea- 
surement will push the phase uncertainty to 3.5/A^. The 
final measurement is then made using the GS state with 
(7 = 1, which lies at the edge of the maximally-squeezed 
Fock regime defined by a <C 1. Thus an arbitrary phase 
can be measured at 3.5/A^ precision with M + 1 measure- 
ments in total. Setting A^m = lO/A/", and solving for M 
gives 



M ^ 0.91n(ln(7Vtan A(9o/2.1)) - 0.4, 



(9) 



where again this is just an estimate subject to run-to- 
run fiuctuations. For Or = tt/S and TV = 10"^, M = 1.6, 
i.e. only 2 or 3 total measurements will be required. 
For N = 10^^, M = 2.6, requiring 3 or 4 measure- 
ments. Even for AOq extremely close to 7r/2, M remains 
small, for example, A^o = 7r/2 — 10~^^ gives M ^ 2.7 
for N = 10^, and M ^ 3.1 for N = 10^^ Hence, 
for arbitrary phases in (— 7r/2,7r/2), our interferometer 
converges quickly to the 3.5/A^ precision within a few 
measurements, regardless of N. The final experimental 
value for the initial unknown phase is then = Xlj=t^ 
with a quantum-limited uncertainty of A^ = 3.5/N. 
The total number of atoms used to obtain this preci- 
sion is Ntot = {M -\- 1)N. For large enough we 
can approximate ln(A^tot/2(M + 1)) ln{Ntot), so that 
M ^1^0.9 \n{\ii{Ntot tan A6>o)) + 0.6, which leads to the 
asymptotic scaling law 



A^ ^ (2.1 + 3.21n(ln(7V,ottanA^o)))/A^tc 



(10) 



That the scaling law should depend on the initial phase 
interval has been previously pointed out [13] . As our final 
approximation effectively overestimates Ntot^ the uncer- 
tainty approaches (p!Q|) from below as N increases. 

In order to verify the accuracy of Eq. ([9]), as well as 
determine the magnitude of the run-to-run fiuctuations, 
we have carried out exact Monte-Carlo simulations of 
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FIG. 3: Monte- Carlo simulation results showing the percent- 
age of runs which achieved the maximum precision of 3.5 /N 
after (M + 1) measurements, plotted versus (M + 1). 



many measurements of the phase 7r/6, with an initial un- 
certainty A^o = 7r/3. During each simulation run, the 
measurement outcome was randomly selected according 
to the output distribution, and the phase information 
was determined numerically via Bayes theorem. The 
prescribed measure-rebalance process was iterated until 
the estimated phase uncertainty reaches 3.5/A^. Figure 
[3] shows the percentage of runs which attain the desired 
resolution on the (M + l)th iteration, for two different N 
values, with 10^ runs each. The averages are M + 1 = 2.2 
for TV = 500, and 2.5 for 5000. Equation ([9]) gives 2.2 and 
2.5 as well. The corresponding variances A(M+1) are 0.4 
and 0.6. We note that as the average approaches a half- 
integer value, the minimum possible variance approaches 
0.5, because M is constrained to integer values. Thus the 
fiuctuations are close to minimum allowed values. 

For TF and PS states, adaptive measurement schemes 
are extremely inefficient. This is due in part to their 
inability to distinguish positive from negative phases, 
which makes rebalancing impossible. But even if this 
were overcome, the primary difficulty is that the phase- 
uncertainty is A^-independent for large phases, so that 
~ N'^ measurements are required to obtain 1/N preci- 

1 /3 

sion. This results in scaling, worse than SQL. 

The above discussions have assumed that the input 
state is optimally squeezed to width amin- A realistic 
input state, however, may deviate from (Jmin^ due to im- 
precise control over u and/or imprecise knowledge of N. 
A straightforward error analysis shows that our scheme 
is extremely robust against such uncertainties. The in- 
crease in the single-measurement phase uncertainty 6 AO 
due to ffuctuations in u and N is found to be 



SAO- 



dAO{N, u) 



dN 



SN ■ 



dh 



5u^ 



(11) 



evaluated at = Umin- The scaling with du'^ reffects the 
fact that u = Umin is a local minimum with respect to 
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the phase uncertainty. This gives 

3M/Mrmn = {2/3)6N/N + {l/8){6u/Urmn)^ (12) 

SO that a 10% variation in N leads to a 7% variation in 
the phase uncertainty, while even a 100% uncertainty in 
u only results in a 13% variation. For our purposes, these 
increases are essentially negligible, and are independent 
of the values of N or 0. Of course there are many other 
potential sources of error, e.g. the precision with which 
the tilt can be rebalanced, and the precision with which 
the scattering length can be set to zero during interfer- 
ometer operation. Reaching the Heisenberg limit in a 
double-well BEC interferometer will clearly require ma- 
jor technological advances in many areas. Assuming that 
a level of precision significantly below the SQL is eventu- 
ally obtained, the scheme we have developed will be the 
optimal method to obtain this precision, whether or not 
it is close to the Heisenberg limit. 

In conclusion, we have shown that an adaptive GS 
state scheme has three advantages over previously dis- 
cussed MZ interferometry schemes. It (1) can readily 
be implemented in a double- well BEC system, (2) can 
achieve a resolution well beyond the SQL for a wide range 
of phases with a single measurement, and (3) quickly 
converges to a final precision ^ 10 /Ntot with only a few 
adaptive measurements. 

This work is supported in part by Nation Science Foun- 
dation Grant No. PHY0653373. 
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